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Abstract 

The effect of Bose-Einstein correlations on multiplicity distribu- 
tions of identical pions is discussed. It is found that these corre- 
lations affect significantly the observed multiplcity distributions, but 
Einstein's condensation is unlikely to be achieved, unless "cold spots", 
i.e. regions, where groups of pions with very small relative momenta 
are produced, occur in high energy heavy-ion collisions. 

1 INTRODUCTION 

Bose-Einstein correlations are important in a variety of physical problems. 
For a recent review cf. In the present report I will review some work on 
such correlations done recently in collaboration with A. Bialas. Let us start 
with some general remarks. 

In particle physics Bose-Einstein correlations imply for identical bosons 
an apparent attraction in momentum space. Since the identical bosons (say 
vr^ mesons) are usually only a subset of the particles present in the final 
state of a high energy collision, it is in general not possible to describe them 
by a wave function. Instead one must use a density matrix p{q,q/), where 
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g = (^1, . . . , qn) denotes the set of momenta of all the identical bosons and 
similarly for q/, or a generalized Wigner function W{P, X) cf. e.g. [0,0,0 
related to the density matrix in the momentum representation by the formula 

p^{P+^,P-^) = J d'^X W{P, X)e-«^. (1) 

If p}^ is the density matrix for N very similar but distinguishable parti- 
cles, the density matrix for the same particles with the additional assumption 
that they are indistinguishable is (cf. 0,0,13] and references given there) 

pf{p.q) = T.PN{P.q.). (2) 

where cr is a permutation of the momenta qi,...,qi\i and the summation 
extends over all the A^! permutations a. Assuming as usual that Trp^^ = 1, 
one finds Trp^ 7^ 1- Thus, the probability of producing identical bosons 
changes from P^{N) to N'P^{N)Trp^, where A/" is a normalization constant. 
This can lead to spectacular changes in the multiplicity distribution 

We limit our discussion to independent particle production models, where 
the multiplicity distribution before taking into account the identity of parti- 
cles is Poissonian, characterized by the parameter u 

P^'\N) = ^e- (3) 

and for a given particle multiplicity A^ the density matrix (for distinguishable 
particles) is a product of identical single particle density matrices 

N 

S\p,'1) = IIPi\p^^^^)- (4) 

i=l 

We make no specific assumptions about the single particle density matrix, 
which is assumed in its most general form 

pf\p,'l) = J2(P\^)^n{n\q). (5) 

n 

Special cases of this model have been considered by many authors. A 
Gaussian pf^ has been particularly popular, because for it the symmetriza- 
tion procedure can be performed explicitly. The original approach P],]^] is. 
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however, quite complicated. We will show that it can be replaced by a very 
simple one, when the problem is translated into the language of thermody- 
namics. 



2 TRANSLATION INTO THERMODYNAM- 
ICS I: DISTINGUISHABLE PARTICLES 

As a warming up exercise and also in order to introduce the notation let us 
consider first the problem of distinguishable particles. The main idea [|10| 



is to interpret the single particle density operator as the canonical density 
operator known from statistical physics 

pf' = T.\n)^n{n\ = \e-^^ (6) 

where the Hamiltonian H = J2n |'^)^n('^| and Z is a normalizing constant. 
By comparison 

An = |e-^- (7) 

and from the condition Trp'i^ = 1 one finds that Z is the standard partition 
function J2n Thus A„ is just the canonical probability that the particle 

is in state n. 

The Hamiltonian corresponding to a given density matrix may be quite 
awkward, but simple cases also occur. E.g. for a Gaussian density matrix 



Pi\p,q) = J-— exp 
V ^vrA^ 



q, K 



2A2 2 



(8) 



where q+ = ^{p + q) and g_ = p — q, the corresponding Hamiltonian is that 
of a harmonic oscillator withQ 



ma; = - (9) 

and 

-'^Thc formula below corresponds to formula (13) of which has been misprinted 
there and should be corrected as give below. 
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1 , 2i?A - 1 

nuj = --log— — — . (10 

There is no problem at RA < \, because from the Heisenberg uncertainty 
principle i?A > ^ and one can check that for i?A ^ also ^ — > oo and hijj 
has a finite limit. Note that according to (10) the frequency of the oscillator 
is temperature dependent. Since, however, we are using thermodynamics 
only as a formalism, this does not cause any trouble. For the same reason 
our use of the canonical distribution does not mean that the system of pions 
is in equilibrium. 



3 TRANSLATION INTO THERMODYNAM- 
ICS II: INDISTINGUISHABLE PARTICLES 

For indistinguishable particles the choice of a single particle as a well-defined 
subsystem is impossible. The procedure recommended by thermodynamics 
textbooks is to choose instead all the particles occupying a single particle 
state. Thus, for each single particle state k there is a subsystem with its 
state defined by the population of state k. The states of the subsystem can 
be labelled by = 0, 1, . . .. The corresponding energies are, of course, n^^fc- 
The subsystem is open, i.e. it can exchange particles with other subsystems. 
Consequently, the probability of state of subsystem k is given by the 
grand-canonical formula 

Pkiuk) = ^e'^(''-'=K^ (11) 

where (i is known as the chemical potential and Zj. - the grand partition 
function - is a normalization constant fixed by the condition Z^^=o Pkin^k) — 
1. By summing the corresponding geometrical progression one finds 

Zk = 1: r. (12) 

The quantity e^^ is known in thermodynamics as fugacity and in our model it 
equals ^. The probability distribution for the whole system can be obtained 
by multiplying the probability distributions for the subsystems 
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P{no,n,,...) = Poino)Pi{m).... (13) 

These formulae yield all that is necessary. We present two examples. 

The average number of particles in subsystem k is given by the textbook 
formula 

, , l91ogZfc uXk ...s 

For vXk small the deviation of the denominator from unity can be neglected, 
and the canonical formula valid for distinguishable particles is reproduced. 
When, however, z/Aq, where Aq denotes the largest eigenvalue of the density 
matrix, tends to one, uq ^ oo, while the average population of the states 
corresponding to A^ < Aq tend to finite limits. This phenomenon is known 
as Einstein's condensation. It was noticed in the Gaussian model by Pratt, 
who called such a system a pion laser 

Let us calculate now the probability that among the pions produced in 
an interaction there are no tt'^-s. The result is of interest, because cosmic ray 
physicists report observations of high multiplicity events without vr'^-s . 
Such events are known as centauros. The probability of no tt'^-s in subsystem 
k is 

Let us notice two limiting cases. For (jik) ^ 1 the formula can be expo- 
nentialized and PkW) ~ e"^"'^^ For irik) 3> 1, -Pfc(O) ^ No 7r°-s in the 
whole system means no vr^-s in every one of the subsystems k. Multiplying 
the corresponding probabilities one finds, when all (jik) <^ 1, 

F(0) ^ e"^, (16) 

where = J2k{f^k)- In the case when almost all the particles are in the k = 

state 

^(0) ^ ^. (17) 

The probability (16) is negligible for of the order of hundred or more, but 
probability (17) is then much larger. It could explain the centauro events, if 
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the TT^-s were condensed. Thus Einstein's condensation could lead to inter- 
esting phenomena (there are more cf. e.g. 0). The problem is, however, are 
we likely to reach condensation in realistic experimental situations. 



4 COLD SPOTS 

Let us start with a simpler question: are the effects of symmetrization on 
the multiplicity distribution large enough to be observable. This question 
must be made quantitative. Suppose the effects are observable, when the 
population of the k = state, corresponding to the largest eigenvalue Aq, is 
doubled due to symmetrization, i.e. when z/Aq > ^. For the Gaussian model 
this condition can be rewritten as 

u>l{RA + ^r. (18) 

For RA = 1 this corresponds to z/ > 1.7. Using formula (14) summed over 
k one finds for 7r+ mesons the equivalent condition A^^+ > 3.4, which for the 
total number of pions yields N,^ ~ 3A^7r+ > 10. This is the number of pions, 
which should be produced in a so-called coherence region, i.e. within a region 
of space-time, where the momenta of the pions produced are sufficiently close 
to each other to make interference effective. As illustrated e.g. by formula 
(8) the contribution from pion pairs, where is large, is negligible. This 
is in fact the reason, why the production regions deduced from studies of 
Bose-Einstein correlations are approximately spherical, while the full pro- 
duction region is believed to be a string, or a bunch of strings. The condition 
Nn > 10 is not very restrictive and, therefore, we do expect the multiplicity 
distribution of pions to be significantly modified by symmetrization. This has 
been confirmed by much more realistic calculations [|12[, [Q, where the un- 



symmetrized distributions had been calculated from LUND models and then 
additional weights corresponding to symmetrization were introduced. This 
observation is interesting, because it shows that the success of the LUND 
model without symmetrization in reproducing multiplicity distributions of 
identical pions is due to a cancellation of errors. First the parameters of the 
unsymmetrized model are adjusted so that it reproduces correctly the ob- 
served multiplicity distribution, which it should not, and then the significant 
correction resulting from symmetrization is not included, so that the final 
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result is satisfactory. 

Let us consider now the question, whether condensation of pions is ex- 
pected to occur in real life. As the quantitative condition we require that 
the k = state contains ten times more particles that the next (triple degen- 
erate according to the Gaussian model) state. In the Gaussian model this 
approximately corresponds to 

iV^+ >30(i?A-^). (19) 

Let us try again RA = 1, which corresponds to Nj^+ > 15. A typical trans- 
verse momentum component for pions is 0.25 GeV. Since A is equal to the 
root-mean-square momentum component, we put A = 0.25 GeV and conse- 
quently i? = 4 GeV"^ ~ 0.8 fm. Thus, in the volume of, say, |7ri?^ ^ 2 fm^ 
there should be on the average more than 15 tt^-s, or more than 45 pions. 
Since the radius of a pion is of the order of one fm, this seems a very large 
density. Alternatively we can estimate the energy density. Each pion carries 



an energy of about ym^ + SA^ ^0.5 GeV. The energy density correspond- 
ing to the formation of the quark-gluon plasma is about IGeV/fm^. Thus 
again we find that the density is too high to be a reasonable starting point 
for a free pion evolution. 

Let us put, however, A = 0.1 GeV. Then R = 2.0 fm and there are 15 
7r+-s per 34 fm^, which is much more reasonable. The corresponding energy 
per pion is about 0.2 GeV. Thus, the energy density is about 0.1 GeV/fm^, 
safely below the density of the plasma phase. If this density is still too high, 
one can go on reducing A and increasing R. We conclude that, if A is as for 
typical pions, the condensation is unlikely. But, if in a heavy ion collision a 
sufficiently large "cold spot" is formed with low A, i.e. with low temperature, 
or equivalently with low relative momenta of the particles produced there, 
there may be condensation of the pions produced in this spot and some 
interesting phenomena may become observable 
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